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We construct the most general effective Lagrangian coupling gravity and electromagnetism up to 
mass dimension 6 by enumerating all possible non-minimal coupling terms respecting both diffeo- 
morphism and gauge invariance. In all, there are only two unique terms after field re-definitions; 
one is known to arise from loop effects in QED while the other is a parity violating term which may 
be generated by weak interactions within the standard model of particle physics. We show that 
neither the cosmological propagation of light nor, contrary to earlier claims, solar system tests of 
General Relativity are useful probes of these terms. These non-minimal couplings of gravity and 
electromagnetism may remain a mystery for the foreseeable future. 



I. INTRODUCTION 

The coupled Maxwell-Einstein system has been exten- 
sively studied. Beyond the classical level, non-minimal 
coupling of the electromagnetic and gravitational inter- 
actions may be described by an effective Lagrangian built 
perturbatively from electromagnetic and geometric ten- 
sors. This is not merely an academic pursuit as such 
terms arc known to be generated within quantum elec- 
trodynamics (QED) by the exchange of virtual charged 
fermions in a curved space-time. 

Specifically, Deser and van Nieuwenhuizen [l| have at- 
tempted to enumerate all possible non-minimal mass di- 
mension 6 actions coupling the Maxwell tensor to the 
Ricmann and Ricci tensors, but excluded parity violat- 
ing ones because QED is a parity conserving theory. 
Following that, Bercnds and Gastmans [2| computed to 
one loop order the photon-photon-graviton 3-point cor- 
relation function, (X '{A^A^ha^}) , for QED in a generic 
weakly curved spacetime; and later, Drummond and 
Hathrell J3(] used their results to do a low energy limit 
"matching calculation" to determine exactly the coeffi- 
cients of the terms obtained by Deser and van Nieuwen- 
huizen. 

With the inclusion of the rest of the electroweak model 
it is conceivable that the weak interactions would induce, 
starting at two loops, parity violating mass dimension 6 
non-minimal terms; although we are not aware of any ex- 
plicit calculation to determine their exact coefficients. It 
cannot be a one loop process because the relevant Fcyn- 
man diagram would need to contain at least one parity 
violating current involving the W or Z boson. 

Because the action is dimensionless' 1 ', any such addi- 
tional terms will be suppressed by some inverse power 
of a mass scale, often times associated with the mass of 
the virtual particlc(s) exchanged. The parity preserving 
terms at mass dimension 6 would receive contributions 



from the standard model beginning at 0(M~ 2 ) and the 
parity violating ones possibly starting at 0(M^ 2 ); where 
M e and M^ are the masses of the electron and W bo- 
son, respectively. These dimension 6 terms could also 
receive contributions from particles that have yet to be 
observed experimentally, if they are massive enough or 
if their interactions arc sufficiently weak. Therefore, one 
may hope that constraining the coefficients of such non- 
minimal terms in the action may in effect probe the exis- 
tence of new physics. Since these are gravitational inter- 
actions, two natural probes are the propagation over cos- 
mological distances of the cosmic microwave background 
(CMB) photons and solar system tests of General Rela- 
tivity (GR). Through explicit calculations, however, we 
will show that one would have to look beyond these av- 
enues to obtain a useful bound. 

In this paper, we shall relax the parity preserving as- 
sumption of Deser and van Nieuwenhuizen and simply 
ask what the entire range of possible couplings between 
the Maxwell tensor and its geometric counterparts is, 
up to mass dimension 6. As we will see, upon field re- 
definition, there arc only two such terms, so that the most 
general form of electrodynamics in curved space-time is 
now given by the action 

S = -^y dWiffl (ft - 2Acc) - \J^xy/^\F^F lw 
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[1] In this paper, the units h = c = fc B = 1 are employed. 



fines the Riemann tensor, and the dual Maxwell tensor 
is defined as F^ = \e il „ a pF al3 . The 1/A X 2 term has a 
well-known contribution from QED. While these effects 
are not new, it is interesting to point out that there is 
really only one non-trivial term. We can compare this 
to the work of Drummond and Hathrell |3j wherein their 
coefficient c/M c 2 = — a EM /(3607rAf c 2 ), where a BM is the 
fine structure constant, is identified with our l/A^. As 
already mentioned, the 1/A 2 2 term may be generated by 



the weak interactions within the standard model. 

It is important to note that the mass dimension 6 ac- 
tions in ((T|) should be viewed as the first terms in an 
infinite series expansion involving the ratios of the micro- 
scopic lengths 1/Ai and I/A2 to either the wavelength of 
the photons or that of the characteristic length scale of 
the gravitation field. That is, the second line of ([lj is 
written down with an implicit assumption that the typi- 
cal energy scale of the photons described by such a the- 
ory has to be significantly lower than Ai and A2. (One 
may see this more explicitly by referring to, for instance, 
the matrix element in equation (2.8) of Drummond and 
Hathrell Q. The mass dimension 6 contributions to the 
energy-momentum tensor Q^ v - i.e. the terms contain- 
ing <?i, <72 and g% - are 0(p 2 /A 2 ) or 0(q 2 /A 2 ) relative to 
the lowest order mass dimension 4 Maxwell contribution 



V, 



fj.i/a/3 



Here p and q are the momenta of the gravita- 



tional field and photon respectively and A 2 ~ M 2 /a EM .) 
Just as real electron-positron pairs could be produced if 
the energy of the photons were of O [few MeV] , one would 
have direct access to the (hypothetical) new physics if the 
energy scale of the photons could reach Ai or A2 and the 
effective theory in ([1} would then no longer be adequate. 
At the very least, one would have to include actions up to 
much higher mass dimensions. As we will see, cosmolog- 
ical and solar system observations are just not sensitive 
enough to constrain a A of the same order of magnitude 
as the photon energies involved (~ 3000K for CMB pho- 
tons and ~ 10 GHz for solar system tests of GR). 

In section |TI] we list the basic tools needed to construct 
all the possible non-minimal terms. Section IlIII includes 
their enumeration from mass dimension two to six, as 
well as an explanation of why, via a re-definition of the 
gauge potential A^ and the metric tensor g^, many of 
these terms are in fact redundant as far as the Maxwell- 
Einstein system is concerned. In section HVl we calculate, 
for cosmological propagation of CMB light and solar sys- 
tem tests of GR, how accurate observations need to be 
in order to place useful bounds on Ai and A2; they ap- 
pear impossible to be achieved in the foreseeable future. 
We conclude and discuss directions for further work in 
section [Vj 



II. BASIC TOOLS 

The most general Lagrangian of gravity and electro- 
magnetism can be written as a sum of the Einstein- 
Hilbert with a cosmological constant (£eh,a cc ) and the 
electromagnetic action (£em) as defined above, and a 
perturbative expansion in the mass dimension of the la- 
grangian density 



/ di x\/ r \g\ f £eh,a cc + £em + jj-C 5 + -j^pC-a + 



with M+ representing the lowest of the (possibly many) 
mass scales that are physically relevant. 



In order to systematically enumerate all the possible 
non-minimal terms coupling gravity to electromagnetism, 
we start by listing the most rudimentary tensors avail- 
able, before forming the general set of scalars out of them. 
The requirements of U(l) gauge invariance and general 
coordinate covariance lead us to the dimensionful tensors, 
the covariant derivative and field strengths 
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with mass dimension 1,2, and 2 respectively. The prim- 
itive dimensionless^ 2 ' geometric objects are 
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where e^ ua P is the fully anti-symmetric Levi-Civita sym- 
bol, and we define e 0123 = — 1. Scalars built from the 
covariant Levi-Civita tensor will violate parity since it 
transforms as a pseudo-tensor (even under parity). As 
far as the complete enumeration of the primitive tensors 
is concerned, the placement of their indices (upper or 
lower) is immaterial because we will be forming scalars 
out of them anyway. The basic strategy for construct- 
ing the most general set of scalars is, then, to consider 
all possible contractions between appropriate products 
of J2]) and ([3]). We need not consider derivatives on the 
e'-tensor as Vt?"" 1 ^ = 0. Furthermore, because 
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we see that it suffices to consider terms that contain zero 
or one extensor only. 

III. ENUMERATION OF TERMS 

Schematically, we seek to form combinations of the 
type 



V a F"R c 



(5) 



As Ffj, v and i? M „ a| a are both of dimension [M] , we 
look for terms satisfying a + 26 + 2c < 6. As all tensors 
except V M have an even number of indices, in order to 
form a scalar we see that a must be an even number, 
therefore no odd mass dimension terms exist. Since we 
arc considering couplings between both the electromag- 
netic and gravitational field, we require b > 0. We have 
summarized these possibilities in Table UJ 



[2] Strictly speaking, the physical dimensions of different compo- 
nents of a given tensor are the same when computed in an 
orthonormal frame, where the metric, in particular, is then 
9nv — > V^v — diag[l, — 1) — lj — !]■ However, since we are forming 
scalars out of these tensors - the result is independent of whether 
the coordinate or orthonormal frame was chosen - we are delin- 
eating this construction within a coordinate frame, where calcu- 
lations are easier. 



Mass 
Dimension 



Possible Combinations of (a,b,c) in ((5J 



(0,1,0) 

(0,1,1), (0,2,0), (2,1,0) 

(0,1,2), (0,2,1), (0,3,0), (2,1,1), (2,2,0), (4,1,0) 



TABLE I: The possible combination of tensors at various di- 
mensions. 



As we enumerate the possible terms, we will not con- 
sider trivial numerical factors as these would be absorbed 
into coefficients of the Lagrangian anyway. The following 
Bianchi identities are quite useful 



v [q f mH = 

Rp.[vap] = 
^[vR\a] P T = 

Using ([7]) it is not hard to show that 
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The anti-symmetric nature of f" implies 

V |1 V„iW' = Rx v F Xv = g^F^ = 
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2. With Levi-Civita Contractions 

Type I: Because of 0, at most 2 indices on the Rie- 
mann tensor may contract with the Levi-Civita, leaving 
only the possibility of F llv R a p€. ,lva ^ , which is zero due to 
the symmetry of the Ricci tensor. 

Type II: The only scalar that can be formed here is 
F^vFapW 11 ^ . Using identity (|6]), we see this is a total 
derivative: 

V M (l^A v F afi ) = F a ?F a p (12) 

Type III: Zero, by the Bianchi identity (|6]). 

B. Dimension: 6 



Here we consider the terms (A) F^Ra^^sRpTcrX, 
(B) F^FprRa^s, (C) F^FapFpr, individual terms 
from (D) \7 p \7 T Fp U R al3l s and (E) \7 p \7 v F al3 F pT . and 
(F) V M VyV p V T F Q/ 3. We again separately consider 
contractions without and with the Levi-Civita tensor. 
As illustrated in the previous section, there are many 
tricks that simplify such constructions greatly. While 
we will not discuss all of them, we have tried to give the 
most salient examples. We have therefore only listed the 
non-zero, non-redundant terms below. 



In 4 dimensional spacctimes, the antisymmctrization of 
more than 4 indices always yields zero. One useful corol- 
lary [4[ is 
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A. Dimension: 2 and 4 

At dimension 2, the only scalar that can be formed here 
is F /JiV g f * u , which is zero. Next we consider the dimension 
4 terms F^Ra^s, F a pF pv , and V^V^F^, which we call 
type 1, II and III, respectively. In addition to contrac- 
tions with the metric, we must separately consider both 
contractions without and with the Levi-Civita tensor. 



Without Levi-Civita Contractions 



1. Without Levi-Civita Contractions 

Type A: All of these are zero because while Fp V is an- 
tisymmetric, all two- index tensors built from contracting 
two Ricmann tensors are symmetric in the indices. As 
an example, using ([7]) it is possible to show that 



Rf T R v 
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where R,^ r R v 



^uppr is (/J <-> v) symmetric. 
Type B: The non- redundant terms are 



f^f^tz 

F a pF^R a ^ 



(14) 



Type I: The only scalar is F^R^ = 0. 
Type II: The scalar this forms is the one found in the 
canonical Maxwell Lagrangian 



V" 



Type III: This is not only a total derivative, which 
does not contribute to the dynamics, it is also identically 
zero from (fTUjl . 



Type C: All of these are zero due to symmetry con- 
siderations. 

Type D: These terms are either equivalent to type A 
(which are zero) or arc of the type VF^y (for some 
appropriate V), which can all be made to vanish at this 
order in mass dimensions by a suitable field re-definition, 
A v -> A v + 8A U , where 8A V ex V v . While it is true that 
this field re-definition will induce a variation of all the 
mass dimension 6 terms, the new terms appearing are of 
even higher order, and so are irrelevant for this analysis. 



Type E: As in [1J, one can show using (|6|) that, up 
to a total derivative, there is an equivalence between the 
term 



(V„*)j 7 ) (V Q F^) 



and 



2 (V Q F Q 7 ) (VpF^) + 2F a J F^ x R Xa + F^F^R, 
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Therefore, the non-zero terms here are equivalent to type 
B and D. 

Type F: This is a total derivative, and does not con- 
tribute to the dynamics. 



2. With Levi-Civita Contractions 

Type A: Zero, for the same reasons as above. To show 
this rigorously, pip may be useful. 

Type B: As an example, using (|TTj) it is possible to 
show 



F aX F 0T R^e^ XT = R^ x F aX F^ - 2R l±a F» x F a x 



(15) 



In all, the non-redundant terms here are 



F^F^K 



F*(3 F pa Ra 



fipa 



(16) 



Type C: Zero or redundant, as above. 

Type D: Zero, as above. 

Type E: These are either redundant or zero, as above. 
Consider the term e^^V^F^VaF^. After integration 
by parts, we find it proportional to 

V a (e^V \F va ) F° p 

= -l^ a F°p(R v aii F Xa + R a a)J F vX ) (17) 

Type F: This is a total derivative, and does not con- 
tribute to the dynamics. 

At this stage, there are six possible non-minimal cou- 
pling terms; three (|14|) of which are consistent with those 
listed by Drummond and Hathrell J3[, while the other 
three (|16|) are new, and are their parity-violating coun- 
terparts. 



C. Absorption of Terms via Metric Re-definition 

Up to this point, the general action is of the form 

S = S'eh,a c 



5 em + I d A x^/\g\ [A^R^ + BK - 



(18) 



where (...) indicates terms that contain neither the Ricci 
tensor nor scalar (but may contain the Riemann tensor). 
We will now demonstrate that A^ v and B can be elim- 
inated through a suitable change of variables, i.e. a re- 
definition of the metric. We make the transformation, 
5g al3 and choose 



g al3 -> g af3 



JW P i 



where A and / Q/3 obey the relations 

Ja[3 = —Aa/3, 



^ 



B 



(19) 



(20) 



where / is the trace of f a p . The overall modification to 
the action flTBl is then 



S = Seh,a oc + "Sem + SS + ... (21) 

SS= f d^/MA^g^ (A MV + B 9lxv ) (22) 

where the (...) represent the same omitted terms. While 
it is true that the change of variables will affect all terms, 
it may be checked that the variation to the non-minimal 
terms occurs at least at mass dimension 8, and so will 
not have an impact on this analysis. 

Discarding the total derivative F^F^, we see that 



SS 



d 4 x 



oi ^ Acc F F» v 

31 M2 V,** 



(23) 



for an appropriate mass scale, M* and constant, /3. This 
rc-definition could pose an issue if it were to flip the over- 
all sign of the Maxwell action, but will not be a problem 
as long as M ^ 3> A cc . This modification to the canonical 
Maxwell Lagrangian can be absorbed by simply rescaling 
the electromagnetic potential as 



A„ 



A n 1 



/3Acc 
AG 



Upon these field re-definitions, we obtain the advertised 
result (P. 



IV. OBSERVABLES AND CONSTRAINTS 

Varying the action ([1]) with respect to A^, the modified 
Maxwell equations arc 

v^ + -i {V[«Rp\u + JWV) F afi 

A l 



^2 (V«^",V/^ + 2F"^ {p R a 



4A 



-2R» vpa V^° ) = 



(24) 



Within the cosmological context, we will derive the 
general solutions to (|2"4"|) in a spatially flat Friedmann- 
Robertson- Walker (FRW) metric and examine their im- 
plications for the propagation of the cosmic microwave 



background. We will then see, as also discussed by J5(, 
although one might have hoped that the large distances 
involved would help accumulate effects from these non- 
minimal terms and render them discernible, the Hubble 
parameter of our universe is simply too small for cosmol- 
ogy to be a sensitive probe. 

We shall also show below that these non-minimal terms 
do modify the geodesies followed by photons in a curved 
background, defining an effective metric, so that the 
travel time and deflection angles of light about massive 
objects will be altered from their standard values. Ob- 
servations of the Shapiro delay of radio signals from the 
satellite Cassini currently provide the tightest bound on 
the PPN parameter 7 PPN . Even though these Cassini 
observations yield tighter restrictions on Ai^ than the 
cosmological ones, they still lie significantly below the 
threshold necessary to provide physically useful con- 
straints. 



A. Cosmological Constraints 

We start first with cosmological probes and work with 
a spatially flat FKW universe, where g M „ = a(ri) 2 r]^. We 
then proceed to solve, in the Coulomb gauge {A$ = 0), 
the general solutions of the vector potential A^ to the 
vacuum modified wave equations, using the JWKB ap- 
proximation. To this end, if k is the spatial momentum 
vector of the photon, it helps to expand the spatial por- 
tion of Ap, A, in terms of basis vectors where one of 
them is parallel to k and the other two correspond to 
left- and right-circular polarizations. That is, if we first 
consider an orthonormal basis defined by unit vectors, 
{ei,en, w\}, then define 



A(r/,x) = (A+ (77) e+ + A_ (77) e_ 



-ikx 



where 



(25) 



(26) 



e± = — = (ei ± ien) 
the resultant equation of motion is (from ([24)) ) 



A'l (1+ V) + A'±ip' + k 2 A ± (1 + x) ± < 2(j)'kA± = 

(27) 

Here k = |k|, the prime denotes derivatives with respect 
to conformal time, 77, and we have used the following 
definitions 



ip(r]) 



A? 



,// 



x(7?) = -a|(^ 



a 1 a 
a 3 \a 2 

>\ 2 



m 



1 l_a^_ 
2A|^3 



(28) 

(29) 
(30) 



Following p, we now attempt a JWKB solution by 
first requiring the solutions take the form 



A±(r]) = „4±exp 



drff±{r}') 



</n 



(31) 



We next assume that time derivatives of A are negligible, 
/' *C f 2 , and proceed to insert (f3i"|) into ([27| . As they are 
small for the cosmological eras of interest, we expand / 
to linear order in <f>, x and 1/) to find (choosing a positive 
root) 

/«fc + *Y + i*(x-VO±0' (32) 



so that 

A± (77) w A exp 



ik(A V +- drl ( X - r(>) (33) 



no 



A-0 

— -±iA6 
2 Y 



The ikAr/ is just the usual plane wave term. The 
parity-conserving 1/A X 2 term contributes to the phase 
a real part, the integral oix~ipj and a dissipativc imagi- 
nary part, iAip/2. The bircfringent ±A0 arises from the 
parity-violating 1/A 2 2 term. 

We may now write 



A± (x exp 



i9 ± - ^AV 



(34) 



with the same proportionality holding for the two cir- 
cular polarizations of the electric field, E±. Hence, 
the energy density of electromagnetic waves propagat- 
ing through the universe will be suppressed by a factor 
e -A0 -jijjg QED contribution can be obtained by bor- 
rowing Drummond and Hathrch's result J3[, that tells 



us that 1/Aj 



-10" 



[/M , which leads us to find 



an extremely small damping of roughly exp [— 10~ 73 ], if 
QED is the most dominant contribution. 

For light coming from a linearly polarized source, over 
the course of its propagation the plane of polarization 
rotates by an angle, |Aa| = \ (6+ —6-) = |A0|. Since 
(in observer time) we have <f>(z) ~ (H(z)/A 2 ) 2 and during 
the matter-dominated era, H 2 ex (z + l) 3 , bounding the 
observed rotation angle restricts the mass scale as 



A 2 >ffn 



(^ + 1) ; 

|Ad 



(35) 



where H w 2 x 10" 33 cV. 

While there are astrophysical (0 < z < 4) sources of 
polarized radiation, such as radio galaxies (see 0] for 
example) which provide 0(1°) limits on polarization ro- 
tation, the CMB turns out to give the best constraint 
because of its large redshift (z m 1100). By rotating 
the plane of linear polarization, birefringence mixes the 
E and B polarization modes of CMB. Specifically, this 



induces a non-zero cross correlation between the tem- 
perature anisotropy and B modes, given by Cj B = 
Cj E sin(2Aa) j|. WMAP @ has put 0(1°) limits on 
the (isotropic) rotation angle of linear polarization of 
the CMB. By applying (|35|) . we obtain a naive con- 
straint of A2 > 10~ 33 McV. However, one must remember 
that the energy of the CMB photons themselves are of 
O(3-3000)K (or, 0(3 x l(T 10 -3 x 1(T 7 ) McV), at least 
23 orders of magnitude greater than this lower bound. To 
obtain a physically meaningful bound, one ought to ask 
instead, how accurate does Aa need to be determined for 
a A2 of at least the same energy scale as that of the pho- 
tons? Setting A 2 > 3000K, and inverting the inequality 
([35]) . the answer is Aa < 0(1(T 55 ). It should be safe to 
assume this is out of observational reach. 

Although this effect is minisculc, it is interesting to find 
a possible standard model source of cosmological birefrin- 
gence that does not invoke any new degrees of freedom 
or extra dynamics (e.g. |10l[ll|). 



To simplify the algebra, it is helpful to rewrite our 
equations in an orthonormal basis using the vierbeins 



/3' 



defined as 



so that 



9^ 



h 



f h vVti 



Key 
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(38) 

(39) 
(40) 



Our conventions here reserve the greek indices for the 
coordinate frame (t, r, #, <j>) and latin for the orthonormal 
frame (i, f , 9, </>). In matrix form, the vierbeins are 



\ 



u 

l/Vu 

r 
\ r sm9 J 



(41) 



B. Solar System Constraints 

Next, we would like to examine how our mass dimen- 
sion 6 terms modify standard GR predictions, so that 
we may use observations to constrain their coefficients. 
In the following, we will employ the JWKB approxima- 
tion to work out the photon's modified dispersion rela- 
tions, due to the addition of the non-minimal terms to 
Maxwell's equations. We then extract the effective metric 
experienced by these photons, and compute the induced 
corrections to both deflection angle and the Shapiro delay 
of light propagating past a massive body. 

Let us first consider the modified Maxwell's equa- 
tions with an example background geometry given by 
the Schwarzschild metric: 

dr 2 = B{r)dt 2 - A(r)dr 2 - r 2 d9 2 - r 2 sin9 2 d(f> 2 (36) 

where B(r) = U(r), A{r) = U~ x {r), with U(r) = 1- ?f, 
r s = 2GM, and M is the mass of the object. Note that 
we have neglected the cosmological constant, A cc . 

Since the wavelength of the light considered here is 
much smaller than the background metric's radius of cur- 
vature, we will use the JWKB ansatz, A^ = Re(a M e l,/i ), 
in which the amplitude a M is slowly varying while the 
phase e 1 ^ varies rapidly. Under these consideration, the 
modified Maxwell equations (|2~4"|) become 

I kftk^dp — -p}R v fj, a pk k a 

(37) 



and e- M is just the inverse of e b u , i.e. 



M„6 



at. 



n»2 ytt-paaP^ "-A"" 



= N v p aP 



For the system of equations to have non-trivial solutions, 
we require det N v a = 0. The eigenvalues of N give us the 
photon's dispersion relations, and the corresponding null 
eigenvectors are the polarization vectors. 



and e a e- M = S a ~. Due to the spherical symmetry of 
the Schwarzschild metric, we can consider, without loss 
of generality, the light propagation to lie in the 6 = | 
plane, i.e. kg = 0. From here on we also choose, for 
simplicity, to assume that the two hypothetical energy 
scales are the same, namely, Ai = A 2 = A. Under these 
considerations, we find 



N b , 



where 



S + A u -A12 -A13 -A 14 

Aia £-A 22 A 23 A 24 

A 13 A 23 £ + A 33 

A w A 24 S + A 44 



(42) 
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r 3 A 2 

(k 2 -k 2 + 2kl)r s 



r3A2 



^22 



M4 



A i2 



r 3 A 2 
kfk^r s 

r 3 A 2 



A 



r 3 A 2 

(k 2 - fc?)r s 

r 3 A 2 



13 



^23 



r 3 A 2 

3^t^<A r s 



r3A 2 



r 



3 A 2 



(43) 



and S = k 2 - kj - k\. 

Two of the eigenvalues of N% are A;? — A;? — k\ = 0; 
this is the canonical light-cone dispersion relation. How- 
ever, their corresponding polarization vectors are pure- 
gauge modes, and thus non-physical. The other two, to 
0(A- 2 ), are 



kf 



ki l ± 



3y/2r a 

A 2 r 3 



(44) 



At this order, if Ai and A2 had been kept distinct, 
the non-minimal modifications to the dispersion relations 
add in quadrature and are thus symmetric under the in- 
terchange of Ai and A2. Therefore, any physical effects 
derived from these relations will not distinguish between 
the two; however, this symmetry does not hold for their 
corresponding polarization vectors. 



which is none other than the geodesic equation. 

Though this effective metric is defined only up to an 
overall conformal factor, if one were to multiply it by any 
function of r, the trajectory of the null geodesies remains 
unaltered. Therefore, the predictions we will quote below 
for the deflection angle and modified Shapiro delay which 
are calculated based on (|45|) are unambiguous. 



Effective metric solution 



2. Deflection Angle 



We wish to analyze the implications of this modified 
dispersion relation in terms of a modified metric, follow- 
ing Myers and Lefrance [12| ] . The two dispersion relations 

(HU) could be viewed as kak~ b g ah = or, in a coordinate 
frame, k^k^g^ 1 " = 0. This defines the effective metric, 
g^v, for each of the two dispersion relations, however, 
only up to an overall conformal factor. We choose to 
write it as 



dr 2 



B{r)dt 2 - A{r)dr 2 - r 2 



where 



B{r)= (1±6) (l 



A(r) = (l±5)[l + - 
\ r 

c 3\/2r s 



A 2 r 3 



(45) 

(46) 
(47) 

(48) 



To 0(r s ), the modification to the dispersion relations 
amount to a re-scaling of B(r) and A(r) found in (|36[) 
by a factor of 1 ± 6, as predicted by the standard result 
in GR. Notice that we continue to work in the 9 = ^ 
plane without loss of generality. Since there is still no 
time dependence, the metric also remains static. 

In this modified geometry the contravariant wave vec- 
tor, fc M = g^k u , is the tangent vector to the path normal 
to the surfaces of constant phase ip, i.e. fc M = dx^ /ds, 
where s is an appropriate affine parameter. In order to 
be convinced that is the right interpretation, all that is 
needed is to show that the modified dispersion relation 
implies that x M (s) is a geodesic of this modified space- 
time. 

We show this, following [13j |. by first taking a covariant 
derivative of the dispersion relation, now with respect to 
the effective metric 



V Q {k^kuW") = 2^V Q fc M = 



(49) 



Since fc M = V^ip = d^ip, it is straightforward to show 
that 



V qKu — V nk a 



Thus 



k V ^k a 



= fc^V^fc" 



(50) 



(51) 



For a general metric of the form (|45|) , the total deflec- 
tion angle of light passing by a massive object is (see, e.g. 



11 ± 

r r 



N 



A(r) 



B(r)r'i 



1 



(52) 



where r$ is taken to be the point of closest approach of 
the light from the object. Expanding the integrand in 
powers of r s /r, the integral above gives us 



A$ 



± 



i\> 



2y/2r 

A 2 r3 



(53) 



The first factor is just the standard contribution from 
the Schwarzschild metric, and the second term is the new 
contribution of the non-minimal terms. An unpolarized 
light ray traversing close to the massive body would incur 
a splitting due to the opposite signs arising in the form of 
(|55|) . While there are observational limits on the deflec- 
tion angle of light passing close to our sun, these limits 
are not as stringent as those derived in the following sec- 
tion, and we will therefore not purse a constraint on our 
non-minimal terms here. The gravitational deflection of 
light computation here suggests that our non-minimal 
terms would also modify the weak lensing signals cur- 
rently sought by large scale structure observations. 



3. Modified Shapiro Time Delay 

We now move on to consider the time-of-flight of a null 
light ray propagating between two points in space such 
that it passes close to a massive object. Such a light ray is 
known to experience a delay in its time of flight, relative 
to the same flight in Euclidean space. This is commonly 
referred to as Shapiro delay [15( . 

For ease of comparison to and discussion in reference to 
the literature, we will now switch to the isotropic gauge 
in calculating the modification to this delay. In this co- 
ordinate system, the Schwarzschild metric (|36p is written 
as 



dr 2 = B(r)dt 2 - A(r) (dr 2 + r 2 d9 2 + r 2 sin 2 



(54) 



where 



B(r) = 
Mr) = 



i _ ik 

_ 4r 

r \ 4 

4r, 



(55) 
(56) 



To first order in r s , both the matrix AT 6 g and the form 
of the modified dispersion relations remains unchanged. 
The expressions that do change are the associated vier- 
beins and the effective metric. In particular, the latter 
becomes 

dr 2 = (1 ± S) B(r)dt 2 - (1 ± S) A(r)dr 2 

-A{r) (r 2 d0 2 + r 2 sin 2 0dcj) 2 ) (57) 



where S was defined in (|46l) . By approximating the null 
path to be a straight line in space, the delay in the round- 
trip time-of-flight between the two points (Pi, 2) is 



At = 2r s log 



4X2X1 
b 2 



± 



iV2r s 
b 2 A 2 



(58) 



If Q is the point on the straight line joining Pi and P2 
closest to the massive object, then b is the distance be- 
tween Q and the object and X\% are the distances from 
Q to Pi. 2, respectively. We note that b is not equal to the 
actual distance of closest approach, r , used in ([52"]) . A 
discussion of this and the different ways of calculating At 
that appear in the literature may be found in appendix 

Within the Parametrized Post-Newtonian (PPN) for- 
mulation, the Schwarzschild metric is altered to quan- 
tify deviations from GR [lfj. To date, the most precise 
measurement of the parameter 7 PPN , which is equal to 1 
in GR, comes from the observation of the Shapiro delay 
from the Cassini spacecraft [17| . Under this parametriza- 
tion, and in the idealized limit where the Earth and 
Cassini are stationary, the Shapiro delay is 



At = (1 +7 P p N )r s log 



4X 2 Xi 
b 2 



(59) 



The Cassini experiment measured the fractional Dopplcr- 
frequency shift of the radio carrier waves, which in turn 
is the time derivative of the Shapiro time delay y(t) = 
dAt/dt. Since the most rapidly changing length scale in 
([59]) is the straight-line closest approach distance, y(t) ex 
db/dt. 

In order to put a constraint on the energy scale A, 
we set y(t) as determined by (|58p equal to y(t) as de- 
termined by ([59]) , thus determining the lower bound on 
A. The actual interpretation and calculation involved 
for the timing measurement is quite involved and the 
reader is referred to [l8( for further information on the 
details on the actual treatment. Using the experimen- 
tal parameter, b ss 6P Q (see both [17| and |19|) and the 
measured value 7 PPN = 1 + (2.1 ± 2.3) x 10~ 5 , we ob- 
tain a naive constraint of A > 0(KU 19 )MeV, which is 



14 orders of magnitude better than the above naive cos- 
mological bound on A2 alone. However, just like in the 
cosmological case, we need to recognize that the energy 
of the radio waves used in the Cassini observations is 
roughly 10 GHz ra 4 x KT 11 MeV, at least 8 orders of 
magnitude greater than the lower bound. Once again, 
we need to ask instead how accurate the timing mea- 
surement needs to be to probe Ai^ ^ 10 GHz. Using the 
second term on the right hand side of ([58]) , the answer 
is At < O(10" 27 )s. This is at least 18 orders of mag- 
nitude more precise than the Casinni observation, if one 
estimates the fractional error for the latter observation 



to be given by the current bound on 7 P 



m 



V. SUMMARY AND DISCUSSION 

We have constructed the most general effective La- 
grangian coupling clcctromagnctism and gravity up to 
mass dimension 6, built from all possible contractions 
between tensors that obey the underlying gauge symme- 
tries of both theories. There arc many such non-minimal 
terms. However, after allowing for field re-definitions of 
the electromagnetic vector potential and the metric, it 
is seen that the number of non-redundant terms reduces 
to two. One represents the type of coupling already ex- 
plored from one- loop quantum effects in QED. The other 
is parity-violating; if it is induced by the standard model, 
we expect it to come from the electroweak sector and to 
be suppressed by 0(M^ V 2 ). 

We have also discussed some of the phenomenology of 
these non-minimal terms, including birefringence of the 
CMB, modified dispersion relations for the photon, as 
well as corrections to the Shapiro time delay, and deflec- 
tion angle. Via detailed calculations, we came to see that 
cosmological and solar system probes do not seem likely, 
within the foreseeable future, to give any physically use- 
ful constraints on 1/Ai and I/A2. This is because, as 
already alluded to in the introduction, observations are 
unlikely to ever reach a level of precision to even probe 
Ai ^-scales of the same magnitude of the photon energies 
involved. 

We end with some suggestions on possible future work. 
Other than the weak lensing surveys already mentioned 
in the body of the paper, Drummond and Hathrell [3j 
have initiated the investigation of these modified photon 
dynamics on a gravitational wave background; it would 



[3] One may consider using the timing measurements of pulsars, 
with masses of order Mq and radii on the order of several km, as 
was exploited by the authors of Q. However, there the measured 
Shapiro delay (really r, the range of Shapiro delay) is governed 
not by the radius of the pulsar itself, as they have claimed, but 
rather the se par ation distance between pulsar and companion 
(see |2(I and 12110 , which is typically of order Rq. Considering 
the relative errors on such observations, solar system observa- 
tions remain a superior test. 



be natural to extend their analysis to include the effects 
of the parity violating 1/A 2 2 term. In this paper, we have 
only examined the dynamics of the photon itself; looking 
at how Einstein's field equations are altered and their 
corresponding implications may provide alternate chan- 
nels to constrain Ai and A2. One may also want to seek 
pcrturbative solutions for A^ or F^ u (with 1/Aj 2 7^ 0) 
about exact solutions of the Einstein-Maxwell system 
(with 1/Aj 2 =0) containing pure magnetic fields, as 
a toy model of more realistic astrophysical systems. Fi- 
nally, a stability analysis of the full system in ([1} may 
also be performed to perhaps help constrain the range of 
physically reasonable values of 1/Aj 2 . 
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Appendix A: Shapiro delay calculations 

For the idealized case where the Earth, satellite and 
Sun are all motionless, the proper Earth-satellite time 
of flight for a light signal (as measured on the Earth) 
has been first computed by Shapiro [15( and elaborated 
in detail in the standard textbook by Weinberg [14[ in 
the standard Schwarzschild gauge. Subsequent calcula- 



tions in isotropic (spatially-conformally-flat) coordinates 
have also been done (see [13], [la], [13 )■ At order r s there 
appear to be differences amongst these calculations, but 
these can be attributed to either the choice of gauge or 
whether or not the straight line approximation is used. 
Reference [19f offers a nice discussion and interpolation 
between some of the methods. 

Most of the calculations found in the literature men- 
tioned above, except in [lj] , use the straight line approx- 
imation to compute the time delay. This method is, in 
fact, exact up to order r S7 at least in the idealized case 
of motionless bodies, because of Fermat's principle in a 
static spacetimc {d t g^ v = got = 0). Namely, the coordi- 
nate time of flight 



At 



gij dx l dxi 
goo ds ds 



ds 



(Al) 



is extremized if x IJ '(s) is a null geodesic of the spacetime 
described by g^ v . In a weakly curved spacetime where r s 
is much smaller than all other length scales, both the null 
geodesies and Ai can be developed as a power series in 
r s . The 0{r s ) accurate Ai can be obtained by employing 
the lowest order solution to the null geodesic equation, 
which is simply a straight line. The contribution to Ai 
due to the deviation of the null path from a straight line 
begins at 0(r%), due to Fermat's principle. 

One could show the equivalence between Weinberg's 
[hi ] and Shapiro's [l5| formulas by an explicit calculation, 
in which the true distance of closest approach, ro, and the 
"straight line" distance of closest approach, b, are related 
through the light deflection angle integral (|52|) . 
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